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We propose a class of qubit networks that admit perfect state transfer of any two-dimensional 
quantum state in a fixed period of time. We further show that such networks can distribute arbitrary 
entangled states between two distant parties, and can, by using such systems in parallel, transmit 
the higher dimensional systems states across the network. Unlike many other schemes for quantum 
computation and communication, these networks do not require qubit couplings to be switched 
on and off. When restricted to 7V-qubit spin networks of identical qubit couplings, we show that 
21og3 A'^ is the maximal perfect communication distance for hypercube geometries. Moreover, if one 
allows fixed but different couplings between the qubits then perfect state transfer can be achieved 
over arbitrarily long distances in a linear chain. This paper expands and extends the work done in 



I. INTRODUCTION 

An important task in quantum information processing 
is the transfer of quantum states from one location (A) 
to another location (B). In a quantum communication 
scenario this is rather explicit, since the goal is the com- 
munication between distant parties A and B (e.g., by 
means of photon transmission). Equally, in the interior 
of quantum computers good communication between dif- 
ferent parts of the system is essential. The need is thus 
to transfer quantum states and generate entanglement 
between different regions contained within the system. 
There are various physical systems which can serve as 
quantum channels, one of them being a quantum spin 
system. This can be generally defined as a collection of 
interacting qubits on a graph, whose dynamics is gov- 
erned by a suitable Hamiltonian, e.g., the Heisenberg or 
Xy-Hamiltonian. One way to accomplish this task is 
by multiple applications of controlled swap operations 
along the communication line. Every external manipu- 
lation, however, inevitably induces noise in the system. 
It is therefore desirable to minimise the amount of exter- 
nal control on the system, to the point that they do not 
require any external manipulation whatsoever. 

Quantum communication over short distances through 
a spin chain, in which adjacent qubits are coupled by 
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equal strength has been studied in detail, and an expres- 
sion for the fidelity of quantum state transfer has been 
obtained HQ. Similarly, in 0, near perfect state trans- 
fer was achieved for uniform couplings provided a spa- 
tially varying magnetic field was introduced. The prop- 
agation of quantum information in rings has also been 
investigated 0. 

In our work we focus on the situation in which state 
transfer is perfect, i.e., the fidelity is unity, and in which 
we can design networks such that this can be achieved 
over arbitrarily long distances. We will also consider 
the case in which no external control is required during 
the state transfer, i.e., we consider the case in which we 
have, after manufacturing the network, no further control 
over its dynamics. In general this will lead us to think 
about more complicated networks than the linear chain 
or chains with pre-engineered nearest-neighbour inter- 
action strengths. We provide two alternative methods 
for understanding how perfect state transfer is achieved 
with pre-engineered couplings. This paper expands and 
extends the work done in Subsequent work has ex- 
amined the extension of this problem to higher excitation 
subspaces The subject of perfect state transfer has 
been independently studied in the first and second exci- 
tation subspaces in where its implementation in an 
array of quantum dots was considered. 

More specifically, we address the problem of arranging 
N interacting qubits in a network such that it allows for 
perfect transfer of any quantum state over the longest 
possible distance. Two qubits are coupled via an XY- 
interaction if an edge connects the two corresponding 
sites. We show further how one can use these networks to 
transfer entangled quantum states and to generate entan- 
glement between distant sites in the network. The con- 
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nection between our approach and the continuous-time 
quantum random walk is highlighted and, in a particu- 
lar example, contrasted with the corresponding result for 
the classical continuous time random walk. 

The paper is organised as follows: In Section ^ we 
set the scene by introducing the problem of perfect state 
transfer in quantum spin systems. In Section IlIII we de- 
rive a necessary condition for our problem in the case of 
graphs with mirror symmetry. This is used in Section Hvl 
to give a limitation of the transfer in chains with uniform 
couplings. In Section we investigate hypercube geom- 
etry as a way to enlarge the previously found limit and 
compare our result to a classical analog in Section IVll 
The quantum walk on the hypercube is found useful to 
derive a modulated spin chain in Section [VIII that allows 
perfect transfer over arbitrary distances. We will exhibit 
a group-theoretic interpretation of this chain in Section 
IVIIII In Sections Ixl and Ixj we consider applications for 
entanglement transfer and the introduction of arbitrary 
phase gates on-the-fly. 



II. STATE TRANSFER IN QUANTUM SPIN 
SYSTEMS 



In order to set the scene, let us first consider quan- 
tum state transfer over a general quantum network. We 
define a general finite quantum network to be a simple, 
connected, finite graph G := {V{G),E{G)}, where V{G) 
denotes the finite set of its vertices and E{G) the set 
of its edges. Two vertices i,j G V{G) are adjacent if 
G E{G). To any such graph G one can associate an 
adjacency matrix A{G) whose elements satisfy 



A,{G) 



lif {i,j)eE{G) 
otherwise. 



(1) 



A quantum spin system associated with such a graph is 
defined by attaching a spin-^ particle to each vertex of 
the graph. To each vertex i g V{G) we can therefore 
associate a Hilbert space Tii ~ C'^. The Hilbert space 
associated with G is then given by 



Hg- (g) 

ieV{G) 



(2) 



where N := 1^^(6)1 denotes the total number of vertices 
in G. 

We define the distance, d(i,j), between any two ver- 
tices V{G) to be the number of edges of the shortest 
path between i and j, i.e., the graph geodesic between the 
two vertices. 

Consider the dynamics of the system to be governed 
by the quantum-mechanical Hamiltonian 



Ha = 



E 



(3) 



We use the symbols erf , erf and erf to denote the familiar 
Pauli matrices acting on the on-site Hilbert space Hi and 
Jij is the coupling strength between the i"* and j*'' sites 
on the graph. Note that Jij = Jji since Hq is Hermitian. 
The total z-component of the spin, given by 



(4) 



(»j)es(G) 



is conserved, i.e., [crf^j, Hg] — 0. Hence the Hilbert space 
Tic decomposes into invariant subspaces, each of which 
is a distinct eigenspace of the operator CTjqj . 

For the purpose of quantum state transfer, it suffices to 
restrict our attention to the iV-dimensional eigenspace of 
crjgj, which corresponds to the eigenvalue (2 — N)/2. Let 
us denote this subspace by Sg- Initial quantum states 
that are in this subspace will remain there under time 
evolution. A basis state in Sg corresponds to a spin 
configuration in which all the spins except one are down 
and one spin is up. Such a basis state can hence be 
denoted by the ket where j is the vertex in G at 
which the spin is up. Thus \j G V{G)} denotes a 
complete set of orthonormal basis vectors spanning Sg- 

When restricted to the subspace Sg, Hg is repre- 
sented, in the above-mentioned basis, by an x iV ma- 
trix which is identical to the adjacency matrix A(G) of 
the underlying graph G H^. The time evolution of the 
system under the action of the Hamiltonian Hg can be in- 
terpreted as a continuous time quantum walk on G (first 
considered by Farhi and Gutmann in 1998 see also 
0). This is because the latter is defined as the time- 
evolution in an A-dimensional Hilbert space spanned by 
states {\j)}, where j G V{G), with a Hamiltonian given 
by the adjacency matrix of G. 

The spin system on the graph G described above plays 
the role of a (noiseless) quantum channel. We see below 
that the continuous time quantum walk on G can be 
viewed as a quantum state transfer along the channel. 

The process of transmitting a quantum state from A 
to B proceeds in four steps: 

1. Initialisation of the spin system to the state |0) :— 
\OaO ■■•00b), which corresponds to the configura- 
tion of all spins down. This state is a zero energy 
eigenstate of Hg [19| . 

2. Creation of the quantum state \'iJj)a& Ha (at ver- 
tex A) which is to be transmitted. Let \tp)A — 
a\0)A+l3\l)A with a,/3 G C and \a\^ + \(3\^ = 1. 

3. Time evolution of the system for an interval of time, 
say to. 

4. Recovery of the state at vertex B, the latter being 
given by the reduced density matrix pg acting on 

The state of the entire spin system after step 2 is given 



3 



by 

|*(< = 0)) = |^^00---00b) 

= a|0^00---OOB)+/?|UOO---OOB) 

= a|0)+/3|l), (5) 

It evolves in time t to 

N 

\^{t))^a\0)+J2f3,{m (6) 

with complex coefficients a,Pj{t), where 1 — jap + 
The initial conditions are given by 
/3^(0) = P and (3j{0) = for all j ^ A. The coeffi- 
cient a does not change in time, as |0) is the zero-energy 
eigenstate of Hg. Hence, it does not even acquire a phase 
factor during the evolution of the state. 

The output state at B after a time t is given by the 
reduced density matrix 




A measure of the overlap between the input state, p^j := 
lip) {ipl, and the output state is the fidehty, 

F{pA,PBit)) : = Tt^pTpbH) p'f 

= y^|a|2 {1-2\PbV+PbP*+P*bP) + \PbV 

(8) 

where it is understood that (3b depends on t. 

Since the |0)^ component of the state is invariant 
under the evolution, it suffices to focus on the evolution of 
the |1)^ component of the state i.e., to the choice a = 
and 13 — 1 in It is therefore convenient to consider 
the transfer fidelity 

fAB{t):^f3B{t)^{B\e-^"^'\A). (9) 

where \A) = |1) = |U0...00i3) and \B) = \N) = 
|OaO . . . 01b) and we have taken h — 1. 

Here we focus only on perfect state transfer. This 
means that we consider the condition 

|/AB(io)| = 1 for some < to < oo (10) 

which we interpret to be the signature of perfect com- 
munication (or perfect state transfer) between A and B 
in time to. The effect of the modulus in IjlOl) is that the 
state at B, after transmission, will no longer be \ip), but 
will be of the form 



The phase factor e^^ is not a problem because (p is inde- 
pendent of a and [3 and will thus be a known quantity for 
the graph, which we can correct for with an appropriate 
phase gate. 

The perfect communication distance d(A, B) is given 
by the distance on the graph, for which perfect state 
transfer is possible. For a fixed number of qubits, N , our 
aim is to find quantum networks which maximise d{A^ B). 
We achieve this in two different ways: 

1. By fixing the nearest-neighbour couplings to be 
identical but considering more complicated graphs 
(see Sections Hvl and IVTl . 

2. By considering linear chains but allowing the 
nearest-neighbour couplings to be different (Sec- 
tions EUl and EIID • 

Note that if there is perfect communication between 
A and i? in a time to, then perfect communication also 
occurs for all times t satisfying 

t^ {2n+ l)to, where n e Z, (12) 

provided the graph is mirror-symmetric (see Section llTHl . 

III. CONDITIONS FOR PERFECT STATE 
TRANSFER IN SYSTEMS WITH MIRROR 
SYMMETRY 

In the rest of the paper, we will examine different 
graphs for the purposes of perfect state transfer. These 
graphs will have mirror symmetry. By mirror symmetry, 
we mean that the graph is identical from the points of 
view of A and B. So, a linear chain with A and B at op- 
posite ends is an example of such a system. The obvious 
question is how can we tell if a proposed graph will per- 
mit perfect state transfer? A necessary condition, as we 
will show, is that the ratios of differences of the eigenval- 
ues of the Hamiltonian, Hq, must be rational provided 
the graph is mirror symmetric. 

With a system capable of perfect state transfer, ini- 
tialised in the state \A), at time to we have the state 

e-'"^*" \A) ^ e"*' \B) (13) 

but by the definition of a symmetric system, A and B 
arc entirely equivalent, and thus after another period of 
time to, we have the state 

and thus the system is periodic, up to a phase 20, with 
period 2to. Thus we conclude that a mirror symmet- 
ric system must be periodic if it is to allow perfect state 
transfer. This may be written most simply as 



a\Q) +e''^(3\l) 



(11) 



I {A\ e-*^G2to \A)\ = \ 



(15) 
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for some time < < oo. 

Let us examine the general state of a periodic system 
with period 2to- We can write 



for eigenstates \ j) of Hq with corresponding eigenvalues 
Ej. Hence for all of the stationary states we have the 
condition 



with fc — 1, 



2E,tQ -2<j) = 2k,n 



(17) 



where the fc^'s are integers. Eliminating (j) between two 
of these, we get that 



(E^ - Ej)2to ^ 2TT{h - kj) 



(18) 



and eliminating the tg between any two of these {E^r ^ 
Ej') gives 



Ei Ej ki kj 



Ei' k/i' ^j' 



(19) 



where Q denotes the set of rational numbers. As the kiS 
are integers, this implies that the ratio is rational. Hence, 
a symmetric system capable of perfect state transfer must 
be periodic, which is equivalent to the requirement that 
the ratios of the differences of the eigenvalues are ratio- 
nal. 



IV. LIMITATIONS FOR PERFECT 
COMMUNICATION OF A UNIFORMLY 
COUPLED CHAIN 

It is desirable to maximise the distance over which 
communication is possible for a fixed number of qubits. 
The optimal arrangement, in this case, is just a linear 
chain of N qubits, where A and B are the qubits at op- 
posite ends of the chain. 

Let us start with the XY chain of qubits, with uni- 
form couplings Ji^i+i = 1 for all 1 < i < — 1. The 
Hamiltonian reads 



H 



(20) 



In this case, one can compute Jab (t) explicitly by diago- 
nalizing the Hamiltonian or the corresponding adjacency 
matrix. The eigenstates and the corresponding eigenval- 
ues are given by 



|fc) 
Ek 



N 



N + 
2 cos 



n=l 

kn 



A + 1' 



nkn 

sm I I in) 



(21) 
(22) 



,N. Thus 



N 



fAB{t) 



N ■ 



Esin ( ^ sin ( 
\N + lJ \ 

k=l 



nkN 
N + 1 



(23) 

Perfect state transfer from one end of the chain to an- 
other is possible for N — 2 and iV = 3, where we find 

2 



-isin(t) and fAsit) 



that fAB{t) 
spectively. 

We have shown that perfect state transfer is possible 
for chains containing 2 or 3 qubits. We will now prove 
that it is not possible to get perfect state transfer for 
iV > 4. 

A chain is symmetric about its centre. Hence the con- 
dition Eq. (|19|l for perfect state transfer applies, i.e. 



Eri 



En 



E , — E 



(24) 



where the Em 's are eigenvalues of the unmodulated chain, 
as given in H22|l . We will explicitly show that there is a 
set of eigenvalues for which this expression does not hold 
for all A^ > 4. 

We are free to choose any values for the indices (pro- 
vided Em' ^ En')^ SO let US choose that m = 2, n = N—1, 
m' = 1 and n' — N . Hence we see, using (|22|) . that we 
require 



2lT 

N+l 



(25) 



N+l 



to hold for perfect state transfer. To find the values of A' 
for which this holds, we make use of algebraic numbers. 
An algebraic number x is a complex number that satisfies 
an equation of the form 



aix 



^ix + an = 0, 



(26) 



with integral coefficients a^. Every algebraic number a 
satisfies a unique polynomial equation of least degree. 
The degree of this polynomial is called the degree of a. 

If a satisfies a monic polynomial (i.e., a polynomial 
with flo = 1) then it is called an algebraic integer of 
degree n. Note that an algebraic integer of degree n is 
also an algebraic number of degree n. Rational numbers 
are algebraic numbers with degree 1, and numbers with 
degree > 2 are irrational. 

Lehmer proved (for example, see T^) that if A^ > 1 
and gcd(/c, A^ -I- 1) = 1, then cos {TTk/{N + 1)) is an al- 
gebraic integer of degree 4'{2[N + l))/2, where is the 
Euler phi-function. 

For n > 3 it can be shown JTjJ that 



e'^ log log n 



(27) 



log log n 



holds, with 7 w 0.5772, Eulcr's constant. Using this 
bound, and by inspection of values not covered by the 
bound, we see that (l){2{N + l))/2 > 3 for A^ > 6. 
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We need to prove that if cos 6 is an algebraic number 
of degree > 3, the quotient 

cos 26* 



cos 9 



is irrational, where 9 — tt/ {N + 1). 

Assume that this expression is rational, i.e. 



cos 29 p 



p,qeZ. 



cosy q 
Using the trigonometric identity 

cos 29 ^2 cos^ 9 - 1 



we can write 



(cos 61)2 - — cos 61 - i = 
^ ^ 2q 2 



(28) 



(29) 



(30) 



which has rational coefficients. According to the def- 
inition, cos 6* is therefore algebraic with degree < 2. 
Given that, from H27|) . cos6' is an algebraic number of 
degree > 3, then we have a contradiction and therefore 
(cos 29) I cos 9 must be irrational. 

Hence we see that for > 6 perfect state transfer 
is impossible because deg(A^) > 3. This simply leaves 
A^ = 4 and N = 5 unproved, which can be done by 
straightforward evaluation. Thus for A^ > 4 

cos 26* , „ 



cos 9 



Hence, perfect state transfer is impossible for unmodu- 
lated chains of length A^ > 4. 



V. PERFECT STATE TRANSFER OVER 
GREATER DISTANCES 

Perfect state transfer over arbitrary distances is impos- 
sible for a simple unmodulated spin chain. Clearly it is 
desirable to find a graph that allows state transfer over 
larger distances, and to that end we examine the d-fold 
Cartesian product of the two-link (three-vertex) chain, 
G. We denote this by C^. 

In general the Cartesian product of two graphs G := 
{V{G),E{G)} andH := {V{H),E{H)} is a graph Gxi/ 
whose vertex set is V{G) x V{H) and two of its vertices 
(g, h) and {g' , h') are adjacent if and only if one of the 
following hold: 

(i) g = g' and {h, h'} e E{H) 

(ii) h = h' and {g,g'} G E{G). 

Let A = (1,1,1,..., 1) and B = (3, 3, 3, . . . , 3) denote 
the antipodal points of G"*. 

We prove that for any dimension d 

\fABit)\ = l, for t^to^^. (31) 

Hence, to is the time for perfect communication between 
the vertices A and B oi G"^. 



Let {X^{G),1 < i < \V{G)\} and {Xj{H),l < j < 
\V{H)\} denote the set of eigenvalues of the graphs G 
and H respectively. The eigenvalues of the adjacency 
matrix of their Cartesian product G x H are precisely 
the numbers: A^(G) + XjiH), with 1 < i < \vIg)\ and 
1 ^ j ^ where each number is obtained as many 

times as its multiplicity as an eigenvalue of the adjacency 
matrix A{G x H). This is because 

A{G xH)^ A{G) ^ lv(H) + ly(G) MH), (32) 

where tv(H) is the |y(i7)| x |V^(iJ)| identity matrix (see 
e.g., 113). 

The eigenvalues of the adjacency matrix of are 
given by 



{jV2|je {0,±l,±2,...,±d}} 



(33) 



and therefore the ratios of differences of the eigenvalues 
are all rational. 

As already observed, the Hamiltonian of a system cou- 
pled via nearest neighbour XY interactions is identical 
to the adjacency matrix. This will hold equally for the 
Cartesian product of G. Hence, 

H^A{G'^) = ^l®^«)y4(G)(g)l^''-^-i (34) 



g-ifft ^ l^-tA(G)t 



(35) 



Thus, if we select a time t — 7r/V2, then we get per- 
fect state transfer along each dimension. Each term 
in the tensor product of applies to a different el- 
ement of the basis (for example, each acts on a differ- 
ent 1 in the definition of ^4, or a different 3 in B). We 
therefore achieve perfect state transfer between A and 
B (as well as between any qubit and its mirror, such as 
(1,1,1,2,3) (3,3,3,2,1)). The fidelity of the state 
transfer is simply the d*"^ power of the fidelity for the 
original chain H23|) . 

Perfect transfer of a single qubit state can also be 
achieved between the antipodes of a one-link hyper- 
cube in any arbitrary dimension, d, in a constant time 
to = 7r/2. This is because perfect transfer occurs across 
a chain of two qubits in this time. 

We can also extend this to the one-link hypercube 
which is coupled via the Heisenberg interaction. This 
is because, in the case of a two-qubit chain, the Hamil- 
tonian in the single-excitation subspace is represented 
by a matrix with identical diagonal elements, and hence 
is the same as the Hamiltonian of an XY model up to 
a constant energy shift, which just adds a global phase 
factor. 

Note that perfect transfer is possible across a ring of 
4 spin-1/2 particles. The topology of this is exactly the 
same as a 2-fold Cartesian product of a one-link chain, 
hence it is a special case of the hypercube we have been 
discussing (whether it is coupled with the Heisenberg or 
XY coupling). 
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VI. CLASSICAL CONTINUOUS TIME 
RANDOM WALK ON THE HYPERCUBE 

In the previous Section we showed that for hypercubes 
generated from both the one-hnk (two-vertex) chain and 
the two-Hnk (three-vertex) chain, the perfect state trans- 
fer time, to, is independent of the dimension d. In this 
Section we will investigate the behaviour of the mean 
hitting time of the classical continuous-time symmetric 
random walk on C^, which we denote by td, and com- 
pare it to io- We will focus our attention on the two-link 
hypercube, since in the quantum case it provides us with 
a greater communication distance, d{A, B) than the one- 
link hypercube does. Unlike tpj we show that td grows 
exponentially with the dimension d (equation (|52(l 'l. We 
also note that the case of the one-link hypercube has 
previously been studied ^3 ■ 

A two-link hypercube is generated by taking the en- 
fold Cartesian product of the graph G := {V{G), E{G)} 
where V{G) = {1,2,3} and E{G) = {{1, 2}, {2, 3}}. 
Hence, the state space for the classical continuous-time 
random walk on is {1, 2, 3}'*. Transitions are allowed 
from a vertex xGG'^to xzte^ where is the i-th unit vec- 
tor and x is a d-dimensional vector with components Xi 
i.e., transitions are allowed to all the nearest-neighbours 
with equal probability. 

If T is the random variable defined as the hitting time 
of B, for a random walk starting at A, then 

tci := IE(r) (36) 

where E(r) denotes the expectation value of T. The 
random variable T can be written as 

N 

T = ^X„ (37) 

i=l 

where iV is a random variable which gives the number of 
jumps that the random walker undergoes in going from 
A to B, and the X^'s are the holding times between suc- 
cessive jumps. We have 

N 

E(T) = E(^X,) 

1=1 

= E{N)E{Xi) 

= E{N), (38) 

where we have made use of the fact that the X^'s are inde- 
pendent and identically distributed with mean E{Xi) ~ 
1. Note that E{N) is the mean hitting time of the cor- 
responding jump chain, which is a discrete-time Markov 
chain. Hence, to estimate the mean hitting time td of 
the continuous-time random walk, it suffices to consider 
the discrete-time random walk given by the jump chain 
of the original walk. 

All the information that we need is contained within 
the 3"^ X 3*^ transition matrix, P . An element Pa,b is 
the probability of transition from A to _B, hence after N 



steps the probability of hitting B (irrespective of whether 
it has previously hit) is (P x P x . . . x P)]^ 3d, which we 

N 

shall denote as P/'gd , and hence 

00 n— 1 

E(A^) = ^ nPi"3d n (39) 

n—2d rn—2d—l 

Since we have a two-link hypercube, it will always take at 
least 2d steps to get from one corner to the opposite one. 
In fact, because we have a two-link hypercube, it will 
always take an even number of steps to hit the opposite 
corner. Thus, P^gd — for odd m and m < 2d. 

So, all we are interested in is the element P^^ 3d , which 
we expect to tend towards a constant value for large n 
i.e., (P^P^")i_3d sa P^^gd as n ^ cx). Thus, we are con- 
fronted with the problem of finding the first element in 
the eigenvector of P^ which has eigenvalue 1 (see, for 
example, [T^'). 

Take, as an example, the special case of c? = 2 where 
the first row of P^" is given by 

( a„ i i i 5„ ) 

The O's occur because in an even number of steps, you 
cannot get to a point that is odd. For n — I, an — 
i (the probability of returning to the start node) and 
bn = (the probability of getting to the exit node). The 
sum of all the elements in the row must be 1 , and hence 
On + bn — ^ for all n. We find a recursion relation for 6„, 
^ + — bn+i, and as n — > 00, we find that 6„ — > i. 

Let us return to the general case. The location of 
a given node of the hypercube can be represented by 
\xi,X2 ■ ■ ■ Xd) where Xi £ {1, 2, 3} i.e., Xi specifies, for the 
dimension, which of the three nodes we are positioned 
on. All of the properties of a given node depend only on 
r and d, where r is a count of the number of Xi = 2. For 
example, the transition rate from one node at r to an 
adjacent node is 

Pncarcst neighbour ~; \ (^0) 

d + r 

This quantity can be understood because the transition 
probability is the same for all connected nodes, and all 
Xi = 1 or = 3 are only connected to Xi = 2 in the i*'' 
dimension, whereas Xi = 2 has two links. 

Given that all properties of a node only depend on 
r and d, that must also be true for the eigenvector, a. 
Hence, we can denote the elements of a by a^. The ele- 
ment of the eigenvector that we are interested in corre- 
sponds to the position |1, 1 . . . 1), i.e., r = 0, so we want 
to find the element a^. We will now find the elements of 
a, the eigenvector of P^. 

In two steps, there are only five ways to get to a specific 
lattice point at a distance r = 2k. 

1. Start at that point, make one step away, and then 
make the same step in reverse. This happens with 
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a probability 



This gives that 



1 



Pretu 



d-2k 



Ak 



d + 2k\d+2k + l d + 2k - 1 



(41) 



2. Start at a distance 2k — 2. There are 4(^2'^) equiva- 
lent jumps, which each occur with probability 

(42) 



{d + 2k~2){d + 2k~ 1)' 



3. Start at a distance 2fc + 2. There are ( 2^ ) '^'^ these 



jumps, each occuring with probability 



(d + 2fc + 2)(rf + 2fc + l)' 



(43) 



4. Start at a distance 2k and go around two edges of 
a square (e.g., |1,2) — > |1,1) — > |2,1)). There are 
4:k{d — 2k) points from which this type of move will 
get us to the specific node we are interested in, and 
this transition happens with a probability 



1 



1 



1 



d + 2k \d + 2k + l d + 2k-l 



(44) 



5. Start at a distance 2k and travel the length of a 
chain (e.g., |1,2) |2,2) |3,2)). There are d- 
2k such jumps, each occuring with probability 



1 



{d+2k){d+2k + l)' 



(45) 



Knowing these, it is possible to write down the ele- 
ments P^a, so we can solve for the eigenvector. 



*V 2 j (d+2fc-2)((i+2fc-l)'^2fe-2 + 

(d-2k\ 2 „ I 

I 2 ) (d+2fe+2)(<i+2fc+l)"2fc+2 + 



(46) 



(4fc+2)(d-2fc) 
d+2k V d+2fe+l 



4fc(d-2fc+l) 

d+2k-l I ^2fc 



a2k- 



Starting with the special case of fc = 0, we see that 

d 



d + 2 



a2. 



and the general solution 



a2k = 



d + 2k 
d + 2k + 2 



0-2k+2 



(47) 



(48) 



can then be proved by induction. We thus have the re- 
quired elements, and just need to normalise them, re- 
membering that there are 2 

0-2k- 



'^^^2k) identical elements 



[d/2] 
fc=0 



-2k 



2k 



a2k 



(49) 



a2k 



d + 2k 
2d3'i-i 



and, in particular 



We know that, as 



ao = f ■ 



(50) 



(51) 



P. 



1,3'' 



gq. Taking this value 
for all n, we can evaluate (|39|) to find that the mean 
hitting time is given by 



tci = E(T) 



00 n — 1 

^ 2nao ]J (1 - ao) 

n=d m—d 

2d-2 + — 



ao 



3^ 



(52) 



The quantum analogue of this mean hitting time is 
given by the time for perfect state transfer between the 
antipodal points A and B. We proved in the previous 
Section that this time is a constant Iq — 7t/\/2. On 
comparing this with (|52|l . we conclude that the graph G"* 
provides an example of a graph for which the quantum 
case leads to an exponential separation. 



VII. PROJECTING A HYPERCUBE ON TO A 
SPIN CHAIN 

Encouraged by the ability of the hypercube to allow 
perfect state transfer (Section 0), we examine the one- 
link hypercube from a different angle. Such a graph falls 
into a general category of graphs, G, that have the prop- 
erty that the vertices can be arranged in columns so that 
there are no edges between the vertices within any col- 
umn, and edges only join vertices in adjacent columns. 
Further, each vertex in column i must have the same 
number of incoming (from column i — 1) and outgoing 
(to column i + 1) edges as all other vertices in that col- 
umn. See Figure ^ for an example. 

Representing the one-link d~dimensional hypercube in 
such a form, we allow the graph G to consist of Nq 
columns. The size of each column (the column occu- 
pation) is given by bi :— \Gi\ — {^^J[^) and the vertices 

in each Gi are labelled Gy , j = {1, . . . , 5J. The i*'' col- 
umn is z — 1 edges away from a corner (say A) of the 
hypercube. 

The only edges are between vertices of adjacent 
columns. From each column there must be a set of edges 
going forwards to the next column, and another set go- 
ing back to the previous one. These are denoted in the 
following manner: 



V!'":^m,,k):je{l,. 



back , 



:{(Gy-,A:) : je{l,. 



,fei},fce{l,. 
,6i},fce{i,. 



.,r,}}(53) 
.,s.}}(54) 
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where and Si denote the number of forward and back- 
ward edges respectively for the i*^ column. Clearly, if all 



the edges are to have ends, \V,^ 



fori 



iq-yback I 
I' i+1 I 



Since there 



is only a single qubit in the first column (61 = 1), each 
vertex in the second column has only a single edge going 
backwards (s2 = 1). With this constraint, and that Si 
and r,: must be integers for all 1 < z < Nq, we require 
that: 



biXi 



bi+lSi+l 

Nc-i 



(55) 
(56) 



The solution that we will choose for this is = Nq — i, 
Si = i — I, which certainly satisfies all conditions. Thus 
we have a graph such that for every pair of numbers (i, j) , 
Gij is connected with Nc — i vertices in Gi+i and each 
vertex in G^+i is connected with i — 1 vertices in Gi. 

Let us define the vectors that span the column space 
■He- 



I col i) 



1 ^' 



(57) 



Farhi et al. 9] note that the evolution with the adja- 
cency matrix Hq of G for this general class of networks 
(not just the hypercube), starting in Gn, always remains 
in the column space Tic because every vertex in column 
i is connected to the same number of vertices in column 
i + 1 and every vertex in column i -f 1 is connected to the 
same number of vertices in column i. 

Thus, we can restrict our attention to the column space 
Tic for the purpose of perfect state transfer from Gn to 
GjVci- The matrix elements of the adjacency matrix of 
G, restricted to this subspace are given by 



J,. 



(col i\ Hg I col i + 1) 

/ Ji 

Ji J2 

J2 J3 

J3 



= ^/i{Nc ~ i) 







Jnc-1 

V Jnc-1 



(58) 



(59) 



This can be seen as follows: 



(col z| Hg |co1 t + l)=^^J2 E^G'^^I \G,+i,r) 



(60) 



Hence, the above graph exhibits the same behaviour as 
the XY chain with "engineered" coupling strengths Ji. 



Nc-l 



(61) 




FIG. 1: An example of a 5-column graph that allows perfect 
state transfer between either end. 



n — 
m = 



1 
-2 
A 



5 
2 
B 



FIG. 2: Couplings J„ that admit perfect state transfer from 
A to i? in a 5-qubit chain. Eigenvalues m of the equivalent 
spin-2 particle are also shown. This is the projection of Figure 
Qon to a chain. 



Such a chain must allow perfect state transfer over any 
length Nc (where |^) = |col 1), \B) = |col Nc)) because 
the hypercube does. In the next Section we prove that 
this is the case using a more physically motivated deriva- 
tion. 

The number of vertices in the graph, G, is given by 
|G| = 2^^'^^, hence it has communication distance of 
log2 |G|. The two-link hypercube in contrast has com- 
munication distance 21og3 |G|. One should note however 
that the degree of each vertex is bounded linearly. 

Some examples of this graph are provided here for dif- 
ferent numbers of columns. 

Nc = 2: two-qubit chain (d=l one-link hypercube) 
Nc — 3: square (d=2 one-link hypercube) 
Nc — 5: for example Figure ^ which reduces to an 
engineered chain, as shown in Figure ||2J). 

For the purpose of perfect state transfer, we have 
stated that the d-dimensional, one-link hypercube is 
equivalent to the graph, G. The equivalence is obvious 
for the case of c? = 1 and d — 2. The general proof arises 
by considering how the Cartesian product of a graph is 
taken when you extend the product from {d — 1) to d 
dimensions. 

Assume the number of vertices in two adjacent columns 
are rii and n^+i in the {d — l)-dimensional hypercube. 
In the {i + 1)*'* column of the d~dimensional hyper- 
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cube, there must still be the Ui+i vertices, plus each 
of the vertices in the previous column have one more 
edge (from taking the Cartesian product). Hence the to- 
tal number of vertices is rii + rii+i. Assuming that the 
{d — 1) -dimensional hypercube has column occupations 
given by a binomial distribution, this specifies that the 
d-dimensional hypercube does as well. Since we know 
that these column occupations hold for d = 2, then by 
induction this must hold for any d. 

What this does not prove is that the edges between 
vertices are correct. This is because they aren't neces- 
sarily correct. While a hypercube must have a specific 
set of edges, the construction of the graph G didn't spec- 
ify which vertices had to be connected to which other 
ones, we just made sure we got the correct number of 
forwards and backwards edges. In that sense, the general 
graph, G, is a 'scrambled' hypercube. No matter what 
this scrambling is, G still reduces to the same chain. 



VIII. STATE TRANSFER OVER ARBITRARY 
DISTANCES 

Suppose we have Nc qubits in a chain, with only one 
qubit in state |t) = |1) and all others in state ||) = |0). 
We previously labelled these as \ denoting that the sin- 
gle excitation is on the j*'* qubit. One may associate a fic- 
titious spin ■^{Nq — 1) particle with this chain and relabel 
the basis vectors as |m), where m — —^{Nc — 1) -I- j — 1, 
as illustrated in Figure |21 This is an equivalent identifi- 
cation to that made in 15], when considering population 
transfer between different atomic levels. 

The input vertex |^) can be labelled as \j = 1) 
or 1 771 = —^{Nc — 1)) and the output vertex \B) as 
\j — Nc) or \m — +^{Nc — I))- Now, consider the 
Hamiltonian, 



appropriate representation of the SU{2) gives 



H = XJx 



Nc-l 

E 



Ji(a, a. 



(62) 



UB(i) = (B\V(i)\A) = 



—I sm 



(64) 



Thus we get perfect transfer of the state from |1) to \Nc) 
in a constant time <o = ""/A. We can select Nc — 1 to be 
divisible by 4 and this eliminates the phase shift caused 
by the factor of —i. 

Note that the case of Nc = 2 is just the same as an 
unmodulated spin chain of the same length, so the calcu- 
lation done previously (|23() is expected to give the same 
result. This it does, provided we remember that in the 
current situation the coupling strength is A/2, whereas it 
was simply set to 1 in the original situation. 

Is there any other inter-qubit interaction in the chain 
that gives Hamiltonian (|59|l when restricted to the single 
excitation subspace? The first choice is the XY model 
with modulated interactions, another one is the Heisen- 
berg model. If we try the Heisenberg model of the form 



(65) 



we obtain 



/ 


^1 


Jl 











Jl 


D2 















■h 


D3 


J3 ... 













J3 


... 















Jnc- 


V 











Jnc~i 


Dnc 



(66) 



where Dj — ^(X]fc=i ^k) — Jj-i ^ Jj- In order to get 
rid of the diagonal elements in the matrix above we can 
apply a magnetic field in the z direction, i.e., we add an 
extra term to H65|) , 



which has the same matrix form as (|59|l . with a scaling 
constant A. 

This corresponds to the flipped spin hopping between 
the vertices j and j + 1 with a probability amplitude 
of Jj. Now, let us choose H to be proportional to the 
angular momentum operator or Jy for some spin J — 
i(iV(7 — 1) particle. In this case the matrix elements 
Jj are i y/j{Nc — j) (these are the same as the elements 
derived in the previous Section up to a numerical factor). 
The evolution of the excitation in the chain is governed 
by the operator 



Nc 



Z^fc=l Jk- 



(67) 



with Bj - \ {Jj-x Jj) - 2(Nc~2) ^k= 

All this means that we can distribute a quantum state 
over any distance with fidelity equal to one as long as 
we engineer the inter-qubit interactions, e.g., the inter- 
qubit distances in the chain, and apply a suitable spa- 
tially varying magnetic field. 



U{t) = exp {—iXt Jx) , 



(63) 



which represents a rotation of the fictitious spin J = 
i(iV(7 — 1) particle. The matrix elements {i'\U\j) are 
well known. Thus working out or looking up an 
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IX. SCALING RELATIONS AND ENERGY 
CONSIDERATIONS 



X. USING THE CHAIN FOR ENTANGLEMENT 
TRANSFER 



In the previous Section we showed that a spin chain 
with engineered interactions can be used to transfer a 
quantum state in fixed time, tg. To compare the com- 
putational complexity of the proposed spin chain, it is 
customary to consider what happens to the energy of the 
system as the number of spins in the chain increases. One 
physical assumption that we might make, for example, is 
that the maximum coupling strength is a fixed size. This 
maximum occurs at the middle of the chain and is 



AJ 



LW/2J 



XN 



Hence, to keep this coupling a constant strength, A must 
scale with 1/iV and to = tt/A must scale with N. 

A second concern is what might happen if we tried to 
extract our state at a time to — St. The fidelity of the 
state transfer is easily approximated from eqn Ht)4|l so for 
small 6t we get: 



fAB{to-St)^l- 



Finally, we could ask the question about what happens 
in the presence of manufacturing errors. In particular, we 
shall consider what happens if the errors only affect the 
eigenvalues of the system. This is not the entire story 
for the spin chain, because we should also consider what 
happens to the eigenvectors (and, in particular, how well 
they maintain their symmetry about the centre of the 
chain since all the eigenvectors are either symmetric or 
antisymmetric). However, in the case of a double appli- 
cation of the chain (which corresponds to nothing hap- 
pening to the stored state), we learned in Section ITTll 
that it is only the eigenvalues that matter. 

Let us assume that we have made some manufacturing 
errors when producing our spin chain i.e., we have some 
errors that are time independent. The ideal energies of 
the eigenstates are Ei and the actual energies are E'^^. 



Jaa = 



{A\ e-'^2to 
\i) 



E 

i 

E 



We can estimate the worst case for the fidelity of the 
identity transformation, by taking the worst error to be 
E'^ — Ei = 6 and by assuming that t^S ^ 1. The error is 
then 



\l-fAA\ 

2to6 



i.e., it scales linearly with N. 



The idea of the rotation of the large spin particle and 
subsequent calculation can also tell us more about the 
system. For example, in the same time that we get per- 
fect state transfer from qubit 1 to Nc, we also get perfect 
state transfer from qubit j to Nc + 1 — j. Under the ac- 
tion of the Jx rotation, these transfers all have the same 
phase. This means that the chain can be used to move 
an entangled state from one end of the chain to another. 
We can start with the BeU state, l/\/2(|01) + |10)), on 
the first two qubits: 



1 



V2 



(|1) + |2)). 



In time to = tt/X this will evolve to the state 



1 



-^{\Nc) + \Nc~l)) 



(68) 



(69) 



having thus transferred the Bell state to the other end of 
the chain. Note that we can not use the state 1 / V2 ( 1 00) -I- 
|11)) because this contains a term with two spins in it, 
and we have restricted ourselves to the subspace of only 
a single spin. We point out, however, that the results 
of show that we will also get perfect state transfer in 
higher excitation subspaces and thus, in principle such a 
state could be transferred. 

The chain can also be used to distribute an entangled 
pair between two distant parties. If we create a Bell state 



1 

71' 



— (|0)^.,|0)c. + |l)Jv/|l>c) 



(70) 



between a non-interacting qubit (NI) and the first qubit 
on the chain (C), then the overall Hamiltonian will be of 
the form 



H' 



H. 



(71) 



Note that the state is exactly the same as the state \i) 
that were were talking about before with the engineered 
chain, but we have to be careful not to confuse those 
states with the states of the non-interacting qubit. The 
state (|70|) then evolves as 



V2 



(|0) 



-iHt 



NI ' 



|0)c 



ii; 



NI 



e-^*ii; 



(72) 



so after the same to, the entangled pair will be the non- 
interacting qubit, and the qubit on the chain. 

1 



V2 



{\0)ni\0)c + ^' 



1 



NI 



(73) 



This prescription is sufficient to transfer the entangle- 
ment of any general two qubit density matrix from being 
between the non-interacting qubit and the f** qubit at 
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t = to being between the non-interacting qubit and the 
Nq' qubit on the chain. This can be understood by see- 
ing how the most general density matrix evolves. What 
we require is that 



TrHG\{A} (p(0)) = Tr-HG\{B} (p(^o)) 
Such a density matrix can be written as 



(ij,i',i')e{o,i} 

Pit) - E 



(74) 



(75) 



a,,,,,,e-''^'* (z'/le'^'*. (76) 



So if a single component of this density matrix evolves, 
giving perfect transfer, so will all the components and 
therefore so will the density matrix as a whole. This 
component evolves as: 



NI 



\3)< 



iH't 



(77) 



INI 



c) 



INI 



((/I 



JHt\ 



After time to, if j or j' were 1, then they will have 
changed to Nc, and if they were 0, they remain as 0. 
Tracing out the effect of all the spins except for the non- 
interacting one and the qubit will return precisely 
the same two qubit density matrix as was initially set up. 
This then allows the density matrix to be split over the 
length of the chain. 

If we want to transmit the complete density matrix, 
we just use two of our engineered chains (Ci and C2) 
in parallel (A^Ci = ^02)- The new Hamiltonian can be 
written as 



(78) 



and an exactly analogous argument now applies so that if 
we create the desired state (which could be the Bell state 
1/V2(|00) + |11)), for example) across the I'** qubits of Ci 
and C2, then after time the state has been perfectly 
transmitted to being on the A^^* qubits of the two chains. 
For an example, see Figure 13 This scheme will work for 
both the engineered spin chain and the hypercubes (since 
the density matrix can be created between the corners of 
two hypercubes). 



Jn - 

C2 



n = 1 
A 



4 
wait 



6 
B 



FIG. 3: Scheme for transferring an arbitrary 2-qubit density 
matrix, p, using two engineered spin chains (Ci and C2). This 
example has a chain length of Nc ~ 6. 



Another choice is to use the Jy rotation (which does 
not give the factor of —i in 1()4|) ). 



Jy — H 



-Ji 
-J2 

-J3 
J3 



'I 








J2 





(79) 



V 







Jnc 



-Jnc- 




Using this in conjunction with the Jx rotation, it is pos- 
sible, along with the transfer of a state through our spin 
chain network, to apply an arbitrary phase gate to it dur- 
ing transmission, simply by choosing the correct linear 
combination of Jx and Jy. Assume that we have picked 
Nc such that Jx gives a phase shift of i. A combination 
of 



jJx ± Vl - j'^Jy 

will thus yield a phase shift e**^ where 

±7 



tan((^) 



(80) 



(81) 



meaning that the initial state {ip) will have evolved to the 
state 



a\0)+e"f'P\Nc). 



(82) 



XI. Jy AND ARBITRARY PHASE GATES 

As previously noted, the Jx rotation introduces a phase 
shift, depending on the length of the chain. There are 
several ways in which this can be avoided. The simplest 
is just to select the correct length of chain. In the case of 
the engineered chain (and also the one-link hypercube), 
if {Nc — I) is divisible by 4, then there is no phase shift 
(since = 1). Similarly with the two-link hypercube, if 
the dimension of the hypercube is even, there is no phase 
shift. 



The final alternative for negating the phase shift, or 
applying an arbitrary phase gate during transmission, 
would be to apply a uniform global magnetic field in the 
^-direction. Applying a field strength B shifts the en- 
ergy of the single spin excitation by B{Nc — 2)/2 and 
the ground state energy is shifted by BNc/'i. Assuming 
transmission of the state occurs in a time t^, then B can 
be selected to give the desired phase shift, (j) by 



B = 



4>~^Nc~l) 
to 



(83) 
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XII. SUMMARY 

We have shown that perfect state transfer is possible 
across a network of qubits, allowing only control over the 
initial design of the network, and no dynamical control. 

When the couplings between adjacent qubits are con- 
strained to be equal, we showed that examples of such 
networks are the one- and two-link d-dimensional hy- 
percube. Perfect state transfer for three- or more-link 
hypercube geometries is shown to be impossible. The 
transfer time is independent of the dimension of the hy- 
percube and for comparative purposes, we calculated the 
expected hitting time in the classical continuous time 
random walk, which increases exponentially with the di- 
mension. 



We have also proposed a spin chain of A'' qubits with 
non-uniform couplings that allows both state and entan- 
glement transfer. This chain can be interpreted in two 
ways: firstly, as a projection of an iV — 1-dimcnsional 
one-link hypercube and secondly, as a rotation in the x- 
direction of a fictitious spin {N — l)/2 particle. 

Finally, we have shown how to effect entanglement 
transfer and how to introduce phases on the transferred 
quantum states on-the-fly. 
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